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I review recent progress in calculating \Vus\ from lattice QCD kaon and hyperon 
systems. A preliminary result from the first dynamical calculation in the hyperon 

o : 

, channel is included. 



I. INTRODUCTION 



^ ' The Standard Model has successfully described strong interactions using a quantum field 

p • theory where universe contains quarks of three generations, which interact via gauge bosons, 
(D ' 

r-| \ the gluons. Although electromagnetism and the strong interaction are not affected by quark 



flavor, the weak interaction may change the flavor of the quarks. This overlap between the 
various generations is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix. 

In 1963, Cabibbo first introduced a 2 x 2 quark mixing matrix to explain semileptonic 
decay in baryons; Kobayashi and Maskawa later extended the matrix to include the then- 



O \ undiscovered bottom quark sector. This becomes the CKM matrix that we are familiar with 
O ■ today: 



^ Vud Vus Vub ^ 



c5 ■ ^CKM 



1) 



Vcd Vcs Vcb 

y Vtd Vts Vtb j 

The Standard Model requires this matrix to be unitary. This gives six unitarity con- 
straints, each of which may be graphically depicted as a unitarity triangle; one such con- 
straint is 

iKdr+iKsr+iK6r = i. (2) 

From the latest PDG 2006jl|, Vud in this unitarity equation is both the dominant term 
and also very well determined from neutron beta decay, 0.97377(27); Vuh is very small. 
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4.21(30) X 10~^. This leaves the middle matrix element Vus as the weak link in deciding 
whether or not this unitarity equation holds. In latest PDG edition, Vus = 0.2257(21) 
is determined to only 0.1%. In this proceeding, I will describe how this number can be 
obtained using calculations from lattice QCD. Note that although currently Vus seems to 
indicate that this unitarity equation holds within error bars, the number has been shifting 
around a lot during the past few years. In 2003, it was 0.2195(23), which is more than 10 
standard deviations away from its current best value. 

In quantum chromodynamics (QCD), physical observables are calculated from the path 
integral. For calculations where the coupling is weak, one can perform the integral by 
hand. However, for long distances the perturbative QCD series no longer converges. Thus, 
to calculate from first principles, one needs help from lattice gauge theory. Lattice QCD 
discretizes space-time such that the path integral over field strengths (especially at strong 
coupling) can be calculated numerically. Since the real world is continuous and infinitely 
large, at the end of the day we will have to take the lattice spacing a — ^ and the volume 
V ^ oo limits to connect to the physical world. However, to simulate at the real pion 
mass (while at the same time keeping the lattice box big enough to avoid finite-volume 
effects) would require much faster supercomputers that have not yet been born. Thus, we 
normally calculate with a few unrealistically large values of the pion mass and then use 
chiral extrapolation to get back to the physical pion mass. 

Here, we quickly mention a few choices of fermion action that have been commonly used 
in lattice QCD calculations. Each has its own pros and cons. They differ primarily by how 
they maintain symmetry, their calculation cost and their discretization error. (Improved) 
Staggered fermions (asqtad)^, ^, 4] are relatively cheap for dynamical fermions, but they 
introduce mixing among parities and flavors or "tastes" , which make baryonic operators a 
nightmare to deal with. The 0(a)-improved Wilson (clover) fermion action is moderate 
in cost and free of the disadvantages of the staggered actions. However, chiral symmetry 
is badly broken at non-zero lattice spacing which causes operator mixing issues. Chiral 

nnriQ 

fermions (e.g. domain-wall (DWF) or overlap) [6|, 17|, l8|, l9[ are free of all the above problems. 
They are automatically 0{a) improved, suitable for spin physics and weak matrix elements. 
These benefits come at great computational cost. Last but not least are mixed actions, where 
one chooses difference valence and sea quark discretizations. Later in this work, we present 
work done using staggered sea quarks (cheap) with domain- wall valence quarks (chiral); we 



FIG. 1: A depiction of the difference between full QCD and the quenched approximation in the 
example of a two-point Green function. 

match the sea Goldstone pion mass to the DWF pion. 

Before we discuss the details of the lattice calculation, it is important to mention the 
"quenched" approximation or 0-fiavor calculations. There are a couple of calculations men- 
tioned later which use this approximation. In the path integral, the correct way to carry 
out the fermionic part of integral is to integrate out the quark/antiquark fields first. This 
leaves the remaining integral as a function only of the gauge links and introduces a fermionic 
determinant. The "quenched" approximation fixes the fermionic determinant as a constant. 
This means that when we calculate, say, a two-point Green function, as depicted in Figure [H 
such as a meson correlator, internal fermion loops have been omitted. 

The quenched approximation was a product of the old days when computers were slow 
and algorithms were not yet sped up to today's standards. Quenching very greatly reduces 
the cost of a lattice calculation by eliminating the fermionic determinant. Of course, modern 
calculations have moved on to focus on urzquenched calculations. This is a lot better, since 
when one ignores the fermion loops, it is very difficult to estimate how much this will affect 
the final numbers; the size of the effect tends to vary greatly between different physical 
quantities. However, since calculations using the quenched approximation can be done very 
fast, one can test new methodologies and ideas using this approximation before performing 
unquenched (dynamical) calculations. The first lattice calculations in and hyperon 
semileptonic decay were demonstrated in this approximation. 

In the following, I will review the latest Vus calculations. There has been a series of works 
devoted to Ki^ decays but only quenched calculations in the hyperon channel so far. In the 
second half of this work, I will show the first lattice dynamical calculation of hyperon decays 
using mixed action. In the final part, I will summarize the current standing of Vus from 
lattice QCD and give some future outlook. 



4 



II. LATTICE Vus CALCULATIONS 

In this work, we will concentrate on three determinations of Vus'- leptonic decay ratios 
(mainly for completion), Ki^^ decays and hyperon decays. We will compare all of them to the 
number listed in PDG 2006. So far, the number from Ki^ decays has the smallest errorbar 
for Vus- 



A. Leptonic Decays 



If one looks at the decays Kf^2 and 7r^2, their branching ratios can be written in terms of 
Vus/Vud and the ratio of the kaon to pion decay constant: 



\Vus 

w 



ud 



Mk{1 



M^(l 



a 



1 + - (Ck - a) 



TT 



ri 



TT 



(3) 



where Ck and C-,^ are the radiative-inclusive electroweak corrections, and the rest of the 
numbers can be obtained from experimental measurements. W. J. Marciano 10|] used decay 
constant ratios jxl ]-k = 1.210(4) (13) from a 2+lf staggered fermion calculation of /tt and fx 
done in 2004 by MILC collaboration and found Vus = 0.2219(25). Of course, there have been 
other fuU-QCD calculations since 2004. For example, RBC/UKQCD use dynamical chiral 
fermions (DWF) to obtain the ratio 1.24(2) ll|. However, none of the other collaborations 



have come out with a number with competitive errorbar yet. In 2006, MILC updated their 



own calculation, 1.208(2)(j4)pj]; this yields = 0.2223 



ill) 



B. Ki3 Decay 

Another way to determine Vus is to look at the Ki^ decay. When one integrates out the 
short- distance dependence, one is left with a low-energy non-perturbative matrix element for 
K to TT, which can be calculated directly in lattice QCD. Using Lorentz invariance, we can 
decompose the matrix element into two form factors (/+ and /„) with differing momentum 
dependence: 

{7r{p')\V,\K{p)) = {p, + p',)U{q') + {p,-p',)f4<l')- (4) 
The so-called "double ratio" technique: 

(7r|g7on|A^)(A^|ll7os|7r) _ , . . 2 xp l"^^ + (r\ 
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FIG. 2: Momentum extrapolation from Ref. 



which has been used to look at B-to-D decays 131], can also be applied to the to vr decay. 
The result is a form factor that only depends on g^, the momentum transfer between the 
initial and final states, and some kinetic factors. This /o can be connected to f± through 

,2 



m 



K 



(6) 



When we extrapolate to = 0, /o = /+. We study different extrapolation forms to estimate 
the systematic error caused by our extrapolation: 



^^2 Quadratic 
,2 \ Polar 



/o(0)(l + Aog^) (7) 

/o(0)(l + Aog^ + cg^) (8) 

= /o(0)/(l-Aog^), (9) 

which are shown in Figure [2] by Becirevic et al. jl^. 

To remove uncertainty due to the momentum extrapolation, we can calculate the matrix 

element directly at = 0. The trick to this technique is to include help from "twisted" 
boundary conditions on the fermions as 



(10) 



The discretized momenta on the lattice will reflect the choice of this 9j as pj = + 
rij^ with rij integer. We can select 9 wisely to cancel out the mass difference, so that 
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we can obtain /+(0) directly. Guadagnoli et al. 15| first demonstrated the advantage in 



the quenched approximation. Later, UKQCD made an exploratory study on a 2+1-flavor 



DWF calculation 161] . They showed that with 50% more statistics, one can get number 
competitive with conventional extrapolation calculations. The advantage of this method 
over the conventional one is smaller or no systematic error due to extrapolation. Thus, 
total error on the calculation is reduced. 

After we obtain /o(0), the next step is to extrapolate the pion mass to the physical one. 



We can get some help from the AdemoUo-Gatto (AG) theorem 171. Il8| . We know that the 
SU (3) symmetry-breaking Hamiltonian is 

H = -^^[rn,-^^yXq. (11) 

The AG theorem tells us that there is no first-order correction due to SU (3)-breaking; thus, 
the correction starts at second order 

/o(0) = /o(0)^^(3) + 0(i72). (12) 

What would be a good measure for SU{?>) breaking? The most natural candidate would 
be the mass splitting between the kaon and pion. So, we expect the remaining correction 
should be small; thus, one would expect the "corrected" lattice /q (after subtracting the 
chiral log), /', should differ from /(f^^'^^ by only a small amount. We construct a ratio 



R{mK,m^) = -jj—^ (13) 



where /"^^(^^ is the S'f/(3)-limit value; in this case, it is 1. We expect the remaining mass 
dependence in Eq. [13] should be relatively small. We then extrapolate the remaining mass 
dependence to the physical sum of the pion and kaon masses-square 

R{mK,m^) =00 + 010^(171^^ + 1711). (14) 

Thus, V„. can then be obtained from 



nKis) = ^|^iK.r^Ewi/+(o)rci4(A.)(i + s^^^,^ + s^^)- (15) 



The decay width, r{Ki3), is taken from experiment, while the phase-space integral /]^, isospin 
breaking (^|[;(2)' long-distance electromagnetic corrections 5^ and short-distance radiative 
corrections Sew are taken from perturbative calculations. 
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Group 




Si 


iV4 (GeV) 


# conf 


/+(0) 


SPQcdR[2^ 





Wilson 


0.500-1.000 


230 


0.961(09) 


JLQCD[22l 


2 


Clover 


0.440-0.960* 


N/A 


0.967(06) 


RBC[23] 


2 


DWF 


0.475-0.700 


94 


0.955(12) 


HPQCD[24] 


2+1 


Staggered 


0.500-0.700 


N/A 


0.962(11) 


RBC/UKQCD[20] 


2+1 


DWF 


0.390-0.700 


150 


0.961(05) 



TABLE I: Summary of existing published /+ calculations from K13 decay 



Table [T] summarizes the results from various lattice QCD groups: quenched, partially 
quenched and full QCD, fermion action variety and the range of the pion mass. Figure [3] 

is taken from individual calculations, combined with the latest PDG 2006 number for 
IZ+V^usI = 0.2169(9). Note the Vus number may be different from the ones given in the 
original papers, due to progress in experimental measurements. The grey band in the graph 
is the range allowed assuming unitarity holds. All the lattice calculations so far have agreed 
with the old estimation from Leutwyler-Roos in 1984[l^. However, not every paper has 
complete estimations of the systematic errors due to lattice artifacts. The work done by 
RBC/UKQCD 20| is one of the exceptions, and thus I quote their number as representative 
of the Vus from Ki^ decay channel: 0.2257(14). 



C. Hyperon Decays 



Hyperon decays provide us with an additional independent channel for determining Vus- 
We start by looking at the low-energy contribution of the transition matrix elements for 
hyperon beta decay, Bi B2e~V; in low-energy effective theory this can be written as 



(16) 



From Lorentz symmetry, we expect the matrix element composed of any two spin- 1/2 nucleon 
states, Bi and B2, to have the general form 



or = /i(g^)7"+ 



M, 



Bi 



M, 



Oi 



9i{q h + — — ^a0q^ + 



M, 



Bi 



Bi 



qa 75 



(17) 
(18) 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

Nf = Wilson 
^ 1 e— 


1 1 1 1 1 


'1 


Becirevic et al. (SPQcdR) 


J. 




Nf = 2 Clover 
' 1 e 1 




Tsutsui et al. (JLQCD) 


Hi 




Nf = 2 DWF 






Dawson et al. (RBC) 


o 


— 1 


Nf = 2+'\ Staggered 






Okamoto et al. (HPQCD) 


— 1 

Hi 




Nf = 2+ I DWF 
^ 1 e- 






Antonio et al. (RBC/UKQCD) 


— 1 




Leutwyler - Roos 






PDG 2006 ' ° 

. 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 


— 1 




,1 


0.22 0.225 




0.23 



FIG. 3: Lattice \Vus\ summary with unified experimental numbers from PDG 

with transfer momentum q = PB2 ^ PBi and V,A indicating the vector and axial currents 
respectively. 

The vector form factor is connected to V^s via 

r = G||K/^(l + 5rad) (19) 



X 



>(20) 



with Am = tub^ — m^j; = Am/mB^-, the radiative corrections 5j 



lad 



and 6g2{fi,gi) taking 

into account the transfer-momentum dependence of /i and gi [25|]. Generally, the ratios 
of gi/fi from experiment and /2//1 in the SU{3) limit are used to get Vus from hyperon 
decays. 



In 2003, Cabibbo et al. 26|] used /2//1 and /i in the SU{3) limit, combined with the 
experimental decay width (or rate) and gi/ fi, to obtain V^g from various channels of hyperon 
decay, as shown in Table [Tll It is not hard to see that if lattice calculations can provide 
better estimates of gi/fi, we can improve the precision of Vus from hyperon decays and 
possibly get a better estimation than the K13 channel. 

So far, there are only two quenched lattice calculations of hyperon beta decay, and they 
are in different channels, S — > n and S*^ S'^. Guadagnoli et al. [^^j extrapolate the matrix 
element S ^ n via an AG ratio, similar to the discussion in the Ki^ decay case, but using a 
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Channel 




(5i//i)"^(^^ 




Tl — 


^ P 


1 n/a 


F + L> 


1.2670(30) 


A - 


Jr 


-^3/2 0.2221(33) 


F + D/3 


0.718(15) 


s- 


— > n 


-1 0.2274(49) 


F - D 


-0.340(17) 




A 


^3/2 0.2367(97) 


F-D/3 


0.25(5) 






072 n/a 


F + D 


n/a 


s° - 


. S+ 


1 0.216(33) 


F + D 


1.32(22) 



TABLE II: Summary of a few hyperon numbers 



dipole form to extrapolate to the zero-transfer momentum point. All of the pion masses are 
larger than 700 MeV; their final numbers are /i = — 0.988(29)stat and Vus = 0.230(5)cxp(7)iat- 



Sasaki et al. 



28 



use lighter pion masses 530-650 MeV and DWF to look at the S° decay 
channel. They extrapolate the vector form factor /i via the variable 6 = {mB2 — ^Bi) /^B2- 
The Ademollo-Gatto theorem suggests the leading-order effect should be 5^, and thus one 
can fit /i(0) to the form cq -|- Ci5'^ . Their final numbers are /i = 0.953(24)stat and Vus = 
0.219(27)exp(5)iat- Unfortunately, the experimental determination of the decay rate is lousy; 
despite /i in H decay channel being compatible within errors, Vus is not well-determined. 
This may further improve in the future with updates from Fermilab KTeV and CERN NA48 
collaborations. One important thing to note is that neither of the calculations has systematic 
error estimates from quenching effects, which we expect might be significant. 

We have taken data looking at both hyperon channels with a dynamical lattice calculation 
for the first time. We use a mixed action, meaning that the sea (staggered) and valence 
(DWF) fermions have different discretization. Our pion masses are relatively lighter than 
the quenched calculations. We only simulate one strange quark mass, which unfortunately 
does not reproduce the correct strange-strange Goldstone mesons. We find a box size around 
2.6 fm and list a few other important parameters in Table IIIII In this work, we report 
on our preliminary calculation in the S~ n channel. We use a projection operator 
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Label 


(MeV) 


niK (MeV) S- 


— > n conf. 


mOlO 


358(2) 


605(2) 


600 


m020 


503(2) 


653(2) 


420 


m030 


599(1) 


688(2) 


561 


m040 


689(2) 


730(2) 


306 



TABLE III: Configuration information 
T = (1 — 7573)(1 + 74)72 in both two- and three-point Green functions and construct a ratio 

^ ^ ^yrg^^c(t.,t,t;,F.,P/; T) / rg|(t,t^,F.; t)~ 

Vrg|(t„t,p,;T)Vr^g(t.,t/,P.;T)' ^ ' 

to cancel out kinetic and overlap Z factors. With multiple insertions of the momentum, we 
can solve for the individual form factors in Eq. [T71 

We need to extrapolate to zero momentum. We use a dipole form, as has been used 
in momentum extrapolation for many baryons' momentum dependence. For the mass ex- 
trapolation, a similar approach to the Ki^ case can be applied here with the help of the 
Ademollo-Gatto theorem. We first construct a ratio and then extrapolate the mass depen- 
dence according to Eq. [131 shown in Figure HI Since our heaviest pion mass is much closer 
to the strange Goldstone meson mass, due to the mass difference in the denominator of the 
ratio, the magnitudes of this point, both central value and especially the errorbar, increase 
a lot. Thus, it provides only a very weak constraint to the fit; this is not an ideal solution 
for our data. 

Alternatively, we can combine the given two-step process into one, by performing a two- 
dimensional fit to the sum and difference of kaon and pion masses, as shown in Figure [51 
But this fit is not ideal either since the constraints from the data points are not strong: only 
four points to define a two-dimensional surface. 

A better constrained fit can be composed by combining the momentum and mass depen- 
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FIG. 4: AG ratio extrapolation to physical + 



dence into a single simultaneous fit: 



1 + (Mi - Mlf (Ai + ^2 (Mi + Ml)) 



(22) 



Mo+Afi(Af2,+M2) 

Figure E] shows the result from simultaneously fitting over all and mass combinations. 
The 2;-direction indicates /i, while the x- and y-axes indicate mass and transfer momentum. 
The surface is the fit using Eq. [22] with color to indicate different masses. The columns 
are the data and the momentum points from different pion masses line up in bands. Our 
preliminary result for fi is —0.88(15). This leads us to a Vus somewhat larger in central 
value than the other calculations but still agrees with them due to the large errorbar. The 
statistics will be greatly improved at the lightest pion mass data in the near future. 



III. CONCLUSION AND OUTLOOK 



To summarize, there are various ways that lattice QCD calculations can help to determine 
Vus in the CKM matrix. (Similar approaches can be applied to the rest of the elements 
with effective lattice fermion actions.) Firstly, we can use the lattice input from the kaon 
and pion decay constant ratios. Currently, MILC has best determined ones, resulting in 
Vus = 0.2226(^5). Secondly, we can use the form factor from Ki^ decay matrix elements: 
here we use the number from RBC/UKQCD, 0.2257(14), in which a sound study and proper 
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FIG. 5: Two-dimensional mass extrapolation after dipole extrapolation to zero-transfer point 

systematics are included. 

Finally, we can use the form factors from hyperon decays. We have started the first 
fuU-QCD 2+1-flavor dynamical calculation. Our preliminary results show consistency with 
previous calculations, but have larger errorbar due to the choice of lighter pion mass. The 
larger statistical error is partially compensated by the decrease in systematic error due to 
extrapolating the pion mass to the physical one. To improve the Vus value from the hyperon 
decays, we need to reduce our statistical error on the vector form factor and improve the 
accuracy on gi/fi to replace the experimental one. Using these strategies, we can make our 
calculation of \Vus\ equivalent to or better than the one from the Ki^ channel. 

Acknowledgements 



HWL thanks collaborator Kostas Orginos for useful discussions on hyperon decays. Com- 
putations were performed using the Chroma software suite 29|] on clusters at Jefferson Lab- 



13 




FIG. 6: Simultaneous extrapolation in and mass 
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